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In this paper some new physical notations are given for the Green's functions and equations of 
motion (EOM) in many body physics with the concept of quasiparticles. It shows how the many 
body correlations existing in many body systems can be taken into account through the derivation of 
EOMs self-consistently. The scheme is compared with the expansion of Feynman diagrams. Finally 
some further discussions are given for the unique properties of EOM techniques in helping to model 
a system and the interactions self-consistently. 
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I. INTRODUCTION 

The concepts and techniques of Green's functions have 
been developed and widely used in the studies of many 
body physics.^ Once the Green's function is known, most 
of the information of the quantum system can be ob- 
tained. Along with the Green's function, the equation of 
motion (EOM) techniques have also been developed. For 
example, different EOM methods have been constructed 
to solve the impurit y pro blems and to solve the many 
body syste m dir ectly^"^ or with the dynamical mean 
field theoryP^EH j n j^s application, the EOM method is 
one of the basic, powerful method to find the Green's 
functions. Moreover, it is also a mono-decisive method to 
find the solution. Under various physical approximations 
it gives a set of closed integral equations which can be 
solved numerically. Therefore, the EOM method is usu- 
ally applied and assumed to be an approximate numeri- 
cal method. However, Comparing with quantum Monte 
Carlo (QMC) method,^ the exact diagonalization (ED) 
method) 17 ! 18 ! and the con tinuo us time quantum Monte 
Carlo (CT-QMC) method, 11212 ^ the EOM method has its 
advantage that it is not only a numerical method but also 
a part of many body theory that can help to develop and 
examine the validity of an existing or a newly emerging 
theory. 

In this paper we will give a few new physical elabora- 
tion of the Green's functions and EOMs, and then discuss 
the relation between EOMs and Feynman diagrams. Fi- 
nally we discuss with example to show how the EOMs 
can verify a Hamiltonian so as to model the system in a 
self-consistent way. 

Now let us imagine an arbitrary correlated many body 
system, as shown in Fig. [T] The particles can be any 
correlated objects. But here in our interest the parti- 
cles will be electrons. All these electrons can be either 
on one impurity site which will then be a single impu- 
rity problem, or on different sites as a lattice problem. 
In a general situation (or presumption) each particle will 
interact with all the other members. If for some cases 
there is no interaction between any pair of members, one 
can still assume that there is interaction between these 
two members but the interaction strength for this 'nonex- 



ist' interaction is zero. Thus all the interactions in the 
whole system can be written as a Hamiltonian %. From 
the many body theory, the complete information of each 
particle can be described with a single particle Green's 
function. However, one should note that here the parti- 
cles are not free particles but with interactions, i.e., they 
are quasiparticles. Therefore, the introduced single parti- 
cle Green's function should be a sum of one term Q 1 ^ naked 
corresponding to the 'naked' particle plus the contribu- 
tions from the quasi-two-body interactions, as illustrated 
m (b) of Fig. [I] We assume that all the three-body in- 
teractions, four-body interactions and other higher or- 
der interactions have been included in this quasi-two- 
body Green's function. Continue this procedure, this 
quasi-two-body Green's function equals a 'naked' two- 
body interaction plus a quasi-three-body Green's func- 
tion as shown in (c) of Fig. [I] One should note that 
there is no isolated two-body interaction or any isolated 
n-body interaction because one n-body-interaction asso- 
ciated members will always interact simultaneous with 
other members outside these n members, i.e., the effect 
of n+1 body interactions and even higher order interac- 
tions will be expressed in this n-body interaction. Repeat 
this procedure again and again, one can clearly observe 
that all the higher interactions/correlations will be cov- 
ered in the quasi-interaction with one order lower, and 
then finally covered by the quasi-single-particle, as pre- 
sented in the following equations, 
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where the indices j,j',j" label difference particles. From 
the above considerations, the achieved single-particle 
Green's functions should be theoretically exact. 

To connect to our theory of EOMs, the above formulae 
just correspond to the derivation of higher and higher 
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FIG. 1: (Color online) Illustration of an arbitrary correlated 
many body system. A) A many body system with interactions 
can be considered as a system of quasiparticles. The wave-like 
lines label the interactions between any two particles. B) A 
quasiparticle can be assumed as sum of a 'naked' particle and 
the relative quasi-two-body interactions, where each ellipse 
circled area labels a quasi-two-body interaction. C) A quasi- 
two-body interaction is similarly a sum of a 'naked' two-body 
interaction and the corresponding quasi-three-body interac- 
tions, where each ellipse circled area labels a quasi-three-body 
interaction. Here in B and C only one quasi-interactions is 
illustrated. Other same order quasi-interactions and those 
with higher order (e.g., quasi- Four-body interactions) can be 
similarly plotted. 

order EOMs as shown in Ref. IMT21 

w « fa\ Si » - ([fa, /]]+>+ « [fa, H];ft »(3) 

W < Ufa' fa] ft > = ([^fa'fa, fa]+) 

+ <[n f a'fa,U];fi^, (4) 



where we have used the double time temperature- 
dependent retarded Green's functions in Zubarev's 
notation,^ /t and f a are the corresponding creation and 
annihilation operators of fcrmions with a and a' label- 
ing different spin channels. Comparing Eqs. ([l]),([2]) with 
([3]) , ([4]) , one can draw a conclusion that the procedure 
of deriving higher and higher order EOMs can distin- 
guish the contributions of the interactions among differ- 
ent number of particles. Therefore, in EOM method it 
is easy to extract the higher order many-body correla- 



tion effects which may be very important in explaining 
the unconventional behaviors of quantum systems. For 
example, if the superconductivity does relate to the pair 
formed with one electron and one hole left by another 
electrons, the higher order long range correlations such 
as fiiafja'fia' must have made great contributions. 

The above EOMs are theoretically exact because each 
Green's function has taken into account all the higher 
order effects. However, when the number of particles in 
this correlated system increases, the system will become 
more and more complicated due to the increasing num- 
ber of EOMs and the number of Green's functions, which 
will make it difficult to explore the system. In some sit- 
uations the EOMs are resolvable, and in some cases they 
are not. For the latter cases, one simple and convenient 
way is to include the most important interactions only. 
Since a higher order interaction will be somehow a minor 
correction to the interaction with one order lower (here 
one should note that the higher order one can never be 
stronger than the interaction with one order lower be- 
cause the interaction with one order lower will cover this 
higher order one), which makes it possible to truncate 
at one order according to the requested accuracy of one 
calculation. If one truncate this procedure at n order, 
a n+ 1 order quasi-particle Green's function will be de- 
coupled, which means that only the n-body interactions 
have been taken into account exactly. All the higher or- 
der interactions beyond n order will be treated with mean 
field approximation which will cause the lose of some fea- 
tures of those correlations with order higher than n. The 
larger the truncation order n is, the more accurate the 
quasi-single-particle Green's function is. At this point 
the EOM method is indeed an approximation method. 
However, one can always get satisfactory precision by in- 
creasing the truncation order though it is a hard work 
to derive the EOMs and made the decoupling. Anyway, 
physical approximations are always used in the studies of 
many body systems, e.g., mean field approximation be- 
tween orbitals, the tight binding approximation between 
sites, or even the approximation with the next or next 
next nearest neighbors. If one increase the truncation 
order of EOMs higher than these physical approxima- 
tions, the numerical error introduced by the decoupling 
can be definitely negligible. 

Next we will compare the EOMs to the diagram expan- 
sion. The usual expansion of the Feynman diagrams is as 
shown in Fig. |2p2l where we only illustrated the Hartree- 
type diagrams. From the expansion, one can observe that 
it is very similar to the EOM method. The higher order 
Feynman diagrams actually just correspond to the con- 
tributions of those higher order correlation effects shown 
in EOMs. However, in physical meaning these two meth- 
ods are quite different. In EOM method the derivation 
of EOMs is actually equivalent to solve the Schrbdinger 
equation (or Dirac equation if relativistically) and in 
each derivation the whole Hamiltonian acts, while dia- 
gram expansion is from perturbation. In EOM method 
each higher order correlation effects will influence and 
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FIG. 2: (Color online) Illustration of diagram expansion. G 
labels a Green's function of a quasiparticle and Go labels 
Green's function of bare particle, a) The expansion of Feyn- 
man diagrams. Here only Hartree-type terms are illustrated 
and complete diagrams should have Fock-type terms, b) Di- 
agrammatic description of Dyson equation, where £ means 
self-energy. 



contribute to the quasi-particle-correlations of one order 
lower, while in diagram expansion it is a simple summa- 
tion of all the orders of diagrams. Therefore, it will lead 
to another difference that each higher Feynman diagram 
can be calculated independent to other diagram, while 
all the Green's functions in EOMs are unknown and can 
only be known by solving the equations self-consistently. 
Moreover, in achieving the Dyson equation, 



(5) 



obviously the higher order terms have also been decou- 
pled in rewriting all the higher order terms as GqTjG. 

From the above discussions, one can learned that, dif- 
ferent to QMC and ED methods, the EOM technique is 
not only a numerical method to solve the eigenvalue prob- 
lem. It can help to explore the many body system, or to 
examine a many body theory. Given a Hamiltonian, if 
the system is well modeled, in the derivation of EOMs all 
the possible quantum many-body correlations will be natu- 
rally reproduced. Otherwise there must be some problems 
or the theory is not self-consistent. For EOM method, 
the Hamiltonian will only consist of the most basic inter- 
actions (or the actual interactions). Those higher order 
correlations whose origin are owing to these basic inter- 
actions will appear in the higher order EOMs but not 
show in Hamiltonian. While for a pure numerical method 
solving eigenvalue problem usually the considered corre- 
lations have to be written in the Hamiltonian, which is 
only one alternative way to study the correlation effects 
because there they can not be self-consistently involved. 
If one want to study more higher correlations, more terms 
will be added into Hamiltonian. In this sense these "nu- 
merically exact" methods can realize and find the solu- 
tion for a give Hamiltonian, and are only numerically 
exact for the given Hamiltonian. Definitely they will not 
exactly reflects the real system if an ill Hamiltonian has 
been given. 



For example, for any a system, the most general model 
should cover all the possible cases and take into account 
all the possible physical effects. Therefore, for a mul- 
tiorbital system the intersite interorbital hoppings and 
on-site interorbital hoppings will be introduced as done 
in Ref . riT)llT2l 

H = — 2^ tijlmfjlafjma + ^ UuhufhiH 
ijlma,i^j il 



^Iraaa' ^ila^ima' 



ilrnaa' ,l<m 

i \^ (y>* ft f ,yl ft f \ (a\ 
~ / < V Ima J ima J ila ~ Ima J ila *> ima ) ' \ J 

ilma,l<m 

The first term on the right hand side (RHS) is the in- 
tersite hoppings. The second and third terms are the 
onsite intraorbital and interorbital Coulomb interaction 
terms. The last two terms are the onsite interorbtial sin- 
gle electron hoppings. The above Hamiltonian is intro- 
duced based on the following presumption. When there 
is no external field and potential, the system should only 
have the kinetic energy and Coulomb interactions, i.e., 
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In the model Hamiltonian the kinetic energy is the sum 
of intersite hopping and onsite interorbital hopping, and 
the Coulomb interactions will exist between any two elec- 
trons so that it should be a sum of intersite and onsite 
Coulomb interactions. However, in usual studies only the 
onsite Coulomb interactions are concentrated, and the in- 
tersite Coulomb interaction are treated as potential with 
mean field approximation and dropped in the Hamilto- 
nian. If directly solving the cluster model or studying 
higher order spacial fluctuations, the intersite Coulomb 
interactions should also be included into the Hamilto- 
nian. 

While in the Kanamori's formj^l the multiorbital 
Hamiltonian is written as 

T~t = — ^ Ujlmfiiyfjma + Ulinu^fHH 
ijlma,i^j il 

"F ^ ^ ^Jlmaa'^ila^ima' 
ilmaa' ,l<m 

— J fii a filer' ftna' fima 

iml.m^l 

^ ] f\l a f\ a i fima' fima i (7) 
iml.m^l 

where the last two terms are the so-called spin-flip term 
and pair-hopping term. Or even some people write the 
Hamiltonian in a more compact form 

7~L = — tijlmfiiafjma 
ijlma,i^j 

+ ^1234 J 'ill fil^fimsfimA- (8) 

i/ml234 
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The spin-flip term and pair-hopping term are usually as- 
sumed to be some kinds of Coulomb interactions and J 
has also the same unit of Coulomb interaction strength 
U. However, the Coulomb interactions can only occur 
between two electrons. In the above equations hf mcr — 
f maimer is the occupation operator of electrons with a 
spin in mth orbital so that Uhf mcr hfi a i means Coulomb 
interaction between two electrons, where I can equal m 
and a can equal a' but can not equal simultaneously. If 
the spin-flip and pair-hopping terms are considered to 
be Coulomb interactions, fj-fmcr must have the physi- 
cal meaning of occupation. However, it labels the onsite 
interorbital single-electron hoppings. Although the ma- 
trix element of f} a fma between two basis functions may 
be nonzero, one should not consider they are Coulomb 
interaction terms. One can also notice that from our in- 
troduction of Eq. (|6| in energy of the system there is also 
no position for these two terms so that they should not 
be added into Hamiltonian. A reasonable way to intro- 
duce the spin-flip term and pair-hopping term is from the 
onsite interorbital single-electron hoppings and consider 
them as some kinds of quantum effect. Apart from the 
onsite single electron hopping. The spin-flip and pair 
hopping terms do relate to the Coulomb interactions. 
The physical meaning of these two terms is that they 
are some kinds of higher order correlation effects between 
two occasional single-electron hoppings, where these two 
hoppings are coupled or correlated with Coulomb inter- 
actions. 

Some people may argue that by choosing special basis 
the onsite interorbital single-electron hoppings may be 
eliminated, which has been discussed in Ref . [T2l but may 
still exist some confusion. Actually this is not always 
true because usually the orthonormal basis are chosen 
to make the matrix element of < 4> f _ l \fi a fm<j\4>u > to be 
zero, while < <fc fJ ,\fi a f m cr\<f > v > is only one of the conse- 
quence of the single-electron hopping term but not the 
interaction itself. Usually the chosen orthonormal basis 
can not guarantee the expectation value of both single- 
electron interorbital hopping term and its two-particle 
offsprings simultaneously to be zero. That's why the 
spin-flip term and pair-hopping can survive in orthonor- 
mal basis. The relation between the onsite interorbital 
single-electron hoppings and the spin-flip term is that 
in EOMs or diagram expansion the spin-flip term can 
be reproduced by the single-electron hoppings. So does 



the pair-hopping term and other two-particle correlations 
shown in Eq. (9) in Ref. [121 If all the onsite interor- 
bital single-electron hoppings vanish (i.e., the interaction 
strengths are all zero due to the physical situation of the 
system, e.g., strictly prohibited by the quantum rules, or 
completely destroyed by the choosing of the basis which 
is a hard or even impossible task when there are many 
partially occupied orbitals), the spin-flip term and pair 
hopping term will also vanish and the Hamiltonian will 
return to 

H = — ^ tijlmfiigfjnuT + Ulfiiltf>>ili 
ijlma,i^j il 

ilmaa' ,l<m 

instead of Hamiltonian in Kanamori's form. On the con- 
trary once spin- flip term and pair- hopping term exist, 
there must be nonzero onsite interorbital single-electron 
hoppings existing, which can be achieved from the EOMs 
in Ref. [T2J When the number of orbital is larger than 
two, these additional two-particle correlation terms in 
Eq. (9) in Ref. H2 have identical amplitudes (or possibil- 
ity) with the spin-flip term and pair hopping term. In a 
self-consistent manner, if one wants to study the effects 
of spin-flip term and pair-hopping term, these additional 
terms in same order can not be neglected. Moreover, one 
may notice that the Hamiltonian in Eq. (|6| can maxi- 
mally have SU(N) symmetry in special physical condi- 
tions where N is the total number of particles in the 
system. As a more general multiorbital model, it can 
of course to describe a special system with only SU(2) 
symmetry by manipulating the prefactors of those in- 
terorbital single electron hopping terms. 

In summary the EOM method is a powerful tool to 
explore a correlated many body system. It can take into 
account the many body correlation effects in deriving the 
higher and higher order EOMs, but not based on the 
perturbation nor adding higher order correlations into 
Hamiltonian. From this point, it can help to build many 
body theory self-consistently and can also have more ad- 
vanced implementations in the studies of correlated many 
body systems. 
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